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ABSTRACT 

We present recent observations from the iJST-Cosmic Origins Spectrograph aimed at character- 
izing the auroral emission from the extrasolar planet HD209458b. We obtained medium-resolution 
(R ~ 20,000) far-ultraviolet (1150 - 1700 A) spectra at both the Phase 0.25 and Phase 0.75 quadra- 
ture positions as well as a stellar baseline measurement at secondary eclipse. This analysis includes a 
catalog of stellar emission lines and a star-subtracted spectrum of the planet. We present an emission 
model for planetary H2 emission, and compare this model to the planetary spectrum. No unam- 
biguously identifiable atomic or molecular features are detected, and upper limits are presented for 
auroral/dayglow line strengths. An orbital velocity cross-correlation analysis finds a statistically sig- 
nificant (3.8 a) feature at +15 (± 20) km s -1 in the rest frame of the planet, at A1582 A. This feature 
is consistent with emission from H 2 B - X (2 - 9) P(4) (\ rest = 1581.11 A), however the physical 
mechanism required to excite this transition is unclear. We compare limits on relative line strengths 
seen in the exoplanet spectrum with models of ultraviolet fluorescence to constrain the atmospheric 
column density of neutral hydrogen between the star and the planetary surface. These results sup- 
port models of short period extrasolar giant planets with weak magnetic fields and extended atomic 
atmospheres. 

Subject headings: stars: planetary systems — stars: individual (HD209458) — stars: atmospheres — 
ultraviolet: stars 



1. INTRODUCTION 

The preceding two decades have witnessed an explo- 
sion in both the number of known extrasolar planets and 
the number of astronomers studying them. There has 
been over an order of magnitude increase in the number 
of detected exoplanets in the last ten years alone, now 
totaling over 400. Among these, transiting planets have 
proven to be the most useful for constraining the physical 
and chemical conditions of the planets and their atmo- 
spheres. HD209458b is a short-period (a = 0.046 AU; 
Henry et a l. 2000) transiting ga s giant, with a mass of 
- 0.7 Mj (|Laughlin et all [2005?) and a density ~ 25% 
that of Jupiter. HD209458b orbits a solar-type (GOV) 
star, which is located at a distance of 47 p c. It was the 
first extrasolar planet observed in transit (jHenrv et al.1 
2000; Ch arbonneau et al"1 120001) . ha s a well determine d 
orbit from radial velocity studies (jMazeh et al.l [2000). 

* Based on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the data archive at 
the Space Telescope Science Institute. STScI is operated by the 
Association of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. 
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and was the first object whose atmosphere was probed by 
absorption line spectroscopy (Charbonneau et al. 2002; 
and see Sing et al. 2008a, b for additional analysis). 

HD209458b was also the first exoplanet whose atmo- 
sphere was investigated usin g far-ultraviolet (far-UV) 
absor ption line spectroscopy (fVidal-Mad iar et al.l 120031 
l200l . The far-UV bandpass is unique in that it offers 
direct access to the strongest transitions of the atoms 
and molecules that constitute the majority of the mass 
in exoplanets (e.g. H, O, C, H 2 , CO, etc.), and that 
are involved in the chemical reactions that produce more 
complex molecules observed in gas giants (H 2 0, CH4, 
C0 2 , TiO, and VO; Swain et al. 2009; Sing et al. 
2009; Desert et al. 2008). Using the f ar-UV capabili- 
ties of HST/STIS, IVidal-Madiar et al.l (12001 detected 
absorption from H I Lya, O I 1304 A, and C II 1335 A. 
These results suggested that HD209458b has an extended 
atmosphere overfilling the Roche lobe and evaporating 
via a hydrodynamic escape mechanism. We note how- 
ever that data analysis techniques and corresponding in- 
terpretations can st rongly influence conclusions drawn 
from these UV data (|Ben- Jaffelll2007l : IVidal-Madiar etHl 
2008). This situation is complicated for O I and C II by 



2 



France et al. 



TABLE 1 
HD209458 COS observing log. 



Dataset 


COS Mode 


Orbital Phase 


T eX p (s) 


Ib4m01 


G160M 


0.00 


4900 


Ib4m02 


G160M 


0.25 


7891 


Ib4m03 


G160M 


0.50 


4912 


Ib4m04 


G160M 


0.75 


7751 


Ib4m05 


G130M 


0.00 


4947 


Ib4m06 


G130M 


0.25 


7942 


Ib4m07 


G130M 


0.50 


4957 


Ib4m08 


G130M 


0.75 


7796 



the low-resolution of the STIS data, and addition al ob- 
servations would be beneficial (jLinskv et al.ll2010j) . 

The majority of studies investigating the atmospheric 
conditions of transiting exoplanets have employed ab- 
sorption techniques. In addition to the UV stud- 
ies cited above, i nfrared (IR) obse rvers have used 
bot h spectroscopic (ISwain et al.l I2009f ) and photomet- 
ric (jD&ert et al.l [2009) methods to constrain the con- 
tent of complex molecules in exoplanet atmospheres. 
Studies of exoplanetary thermal emission in the mid- 
IR have shown that these objects can be observed di- 
rectly at wavelengt hs where the re l ative planet /star con- 
trast is favorable (jKnutson et al.l 120081 I2009D . Analo- 
gously, both Jupiter and Saturn are strong far-UV emit- 
ters, well above any contribution from scattered sunlight 
at A < 1700 A (|Feldman et al.|[l99l[Gustin et al.l [20021 
120091) . 

Far-UV emission from Jupiter and Saturn is domi- 
nated by atomic (Lyman series) and molecular hydrogen 
(H2) lines, the primary c onstituents of gas giant plan- 
ets (jSudarskv et al.l 120031) . H 2 may be responsible for 
the Rayleigh scattering obse rved in optical spectra of the 
HD20 9458b atmosphere (jLecavelier Pes Etangs et al.l 
2008). H2 in the atmosphere of giant planets is excited 
by two primary mechanisms: electron bombardment and 
fluorescence pumped by stellar emission lines. Both pro- 
cesses excite the ambient molecules to higher-lying elec- 
tronic states whose decay produces a highly structured 
spectrum in the 700 - 1650 A bandpass. The majority of 
these lines can be attributed to fluorescence from the Ly- 
man and Werner bands of H 2 (-B x £+ andC 1 ^ toX^j, 
however higher electronic states contribute to the ob- 
served electron impact spectrum). Gas giant aurorae 
produ ce the highest surface brightnesses (jClarke et al.l 
1998) as the electron impact spectrum is most intense 
where the magnetic field lines connect with the atmo- 
sphere. Dayglow (non-auroral illuminated disk) emis- 
sion from the giant planets also shows strong features of 
H2 , with differing levels of molecul ar excitation and con- 
tribution from solar fluo rescence (jFeldman et al.l 119931 : 
IWolven fc Feldmanl fl998l) . It follows that the far-UV 
could be an ideal wavelength regime in which to probe 
the atmo spheres of hot Jupiters through their direct 
emission (Yclle 2004). Most transiting extrasolar giant 
planets are at relatively small separations from their par- 
ent star (a < 0.1 AU), where the line flux from the host 
star should be more intense than at Jupiter by a large 
factor (> 10 3 ) and bombardment from a stellar wind 
may be enhanced relative to the Jovian environment. 
Additionally, by virtue of their tight orbits, these plan- 



ets have large orbital velocities (> 100 km s _1 ). Hence, 
by observing the system near both quadrature positions 
(Phase 0.25 and Phase 0.75), even with modest spectral 
resolution, the observed velocity shift between the two 
epochs can provide confirmation that the signal is plane- 
tary in origin. With a velocity resolution of ~ 15 km s _1 
and an order of magnitude increase in sensitivity over 
previous far-UV instruments, the Cosmic Origins Spec- 
trograph (COS) enables new opportunities for studies of 
exoplanetary atmospheres. 

In this paper, we present far-UV observations designed 
to search for direct emission from an extrasolar planet. 
We use early observations from COS, recently installed 
on HST, to constrain the emission properties of bulk 
atmospheric constituents such as H2, as well as several 
atomic species. We describe the COS observations and 
data reduction in Section 2. In Section 3, we describe the 
details of an exoplanet auroral emission model based on 
observations of Jupiter, and present a quantitative anal- 
ysis of the far-UV spectrum of HD209458(b) obtained 
with COS. In Section 4, we compare the data to syn- 
thetic spectra, and place limits on the H2 emission in 
the HD209458b atmosphere. We conclude with a brief 
summary in Section 5. 

2. HST-COS OBSERVATIONS AND DATA REDUCTION 

The Cosmic Origins Spectrograph is a fourth- 
generation instrument installed on HST in 2009 May 
as part of Servicing Mission 4, STS-125. The spec- 
trograph is optimized for observations of point sources 
in the far-UV bandpass (1150 < A < 1750 A) and 
achieves high sensitivity in the far-UV by using a sin- 
gle optical element for dispersion, focus, and correction 
of the spherical aberration of the Hubble primary mir- 
ror. COS covers the 1150 - 1750 A band at moder- 
ate resolving power (G130M and G160M; R ~ 20,000; 
Av « 15 km s _1 for point source observations), with 
additional modes that offer low resolution spectroscopy 
to wavelengths as short as the Lyman Limit (G140L; 
R ~ 2000; 900 < A < 2100 A; McCandliss et al. 2009) 
as well as low and medium-resolution spectroscopy in 
the near-UV (G185M, G225M, G285M, and G230L; 
1700 < A < 3200 A) First COS science results can b e 
found in lFrance et~atl (|2009) and lDanforth et al.l (j2009D , 
and a full instrument description and on-orbit perfor- 
mance characteristics are in preparation (J. C. Green et 
al. 2010, in preparation; S. Osterman et al. 2010, in 
preparation) . 

HD209458 was the first science target to be ob- 
served with COS after nominal science operations be- 
gan in 2009 September (program 11534, HST dataset 
lb4m). The observing program was divided into eight 
visits, four phases of the exoplanet orbit in G130M and 
G160M. The relevant observational parameters are given 
in Table 1, and we used the o rbital period and zero 
phase data from lKnutson et al.l (120071 ) . P = 3.52474859 
days, T c = 2,452,826.628521, to determine the orbital 
phase during the observations. The transit duration for 
HD209458b is ~ 3.2 hrs, thus a single transit spans 
approximately two HST orbits. The transit data will 
be presented in Linsky et al. (2010), and we refer the 
reader to this work for a discussion of atomic absorp- 
tion from the extended atmosphere of HD209458b. In 
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Fig. 1. — The far-UV spectrum of HD209458. This plot is a weighted coaddition of Phases 0.25, 0.5, and 0.75, combining data from the 
G130M and G160M modes of HST/COS. The flux units are those preferred in the COS instrument handbook, the femto-erg flux unit, 
FEFU (1 FEFU = 1 x 10~ 15 ergs cm" 2 s" 1 A" 1 ). 



the present work, we will focus on the observations made 
at the Phase 0.25 and Phase 0.75 quadrature positions, 
using the Phase 0.5 secondary eclipse of the planet by 
the star (hereafter "occultation" ) observations as a base- 
line reference for any emission signal attributable to the 
planet. We present the 1150 - 1650 A spectrum in Figure 
1. 

HD209458 was observed by HST-COS in the interval 
from 15 September 2009 to 18 October 2009, for a to- 
tal of 20 orbits. All observations were centered on the 
host star (R.A. = 22 h 03 m 10.82 s , Dec. = +18°53'03.8"; 
J2000) and COS performed target acquisition using near- 
UV imaging through the Bright Object Aperture (BOA). 
In order to achieve continuous spectral coverage and min- 
imize fixed pattern noise, observations in each grating 
were made at four central wavelength settings (A1291, 
1300, 1309, and 1318 in G130M and A1589, 1600, 1611, 
and 1623 in G160M). The data were processed with the 
COS calibration pipeline, CALCOS0, v2.11. During the 
analysis, we observed anomalous behavior in several stel- 
lar emission lines that was traced back to a change in flux 
calibration reference files over the course of the observ- 
ing period. The entire data set was reprocessed following 
completion of the program using the most recenlQ COS 
sensitivity files. 

3. HD209458B AURORAL MODEL AND SPECTRAL 
ANALYSIS 

3.1. Atmospheric H% Emission Model 

In this section, we present an H2 model for the atmo- 
spheres of giant planets. We describe the model in gen- 

3 We refer the reader to the cycle 18 COS Instrument 
Handbook for more details: http://www.stsci.edu/hst/cos/ 
documents/handbooks/ current/cos_cover .html 
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eral terms as a tool for future far-UV analyses of exoplan- 
etary dayglow/aurorae, and use this model to compare 
with the COS observations in Section 4. Synthetic spec- 
tra of fluorescent emission from H2 can be made by com- 
puting the collisional and radiative excitation rates into 
the upper electronic states. This model assumes a ground 
electronic state population, then uses excitation cross 
sections and an incident radiation field to calculate the 
population distribution in the rovibrational levels of the 
upper electronic state (predominantly B l Y,+ and C 1 ]!^). 
The molecules will then return to the ground electronic 
state following the appropriate branching ratios, produc- 
ing far-UV emission lines (700 < A < 1650 A) and leav- 
ing the molecules in excited rovibrational levels. We have 
adopted a modified version of the sy nthetic molecu l ar hy - 
drogen emission code presented by IWolven et al.l ()1997f ) 
to model fluorescence induced by electron impact and so- 
lar Lya at the Shoemaker-Levy 9 impact site on Jupiter 
(see also models of the far-UV Jovian dayglow presented 
by Liu & Dalgarno 1996). 

We begin by creating an electronic ground state pop- 
ulation for the molecules. This is done following the for- 
malism of Herzberg to create a matrix of possible en- 
ergy levels, characterized by their rotational and vibra- 
tiona l quantum numbers, J and v, respectively (jHerzbergl 
1950). An occupation probability of a given rovibrational 
state, P{v, J), is found multiplying the Boltzmann factor 
by the degeneracy of a given state. 

The photon-induced fluorescence starts from the 
ground state population, the molecules are then "illu- 
minated" by the EUV97 model solar radiation field (Fq 
from 1 - 2000 A; Tobiska & Eparvier 1998). Excitation 
into the upper states is determined by theoretical rates 
at the Born limit, using the transition probabilities of 
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Abgrall et al. (1993a, b). The excitation rate is given by 



R(vJ^nv'J') = h 



(1) 



where v and J are the vibrational and rotational quan- 
tum numbers of the ground electronic state and v' and J' 
are the vibrational and rotational quantum numbers of 
the electronic level n. The units of R and I\ are photons 
cm~ 2 s _1 A -1 . r\ is the total optical depth at a given 
wavelength, and r™ is the optical depth for the transi- 
tion connecting vJ and nv'J'. The intensity of radiation 
being absorbed at a given wavelength, I\, is 



h 







(2) 



where 1° is the model solar spectrum, with a scaling 
proportional to the HD209458 - HD209458b distance. 

Following the absorption, the molecules emit obeying 
branching ratios. The emission rate from a particular 
level of nv'J' to a level v" J" of the ground electronic 
state is given by 



j~) ( / t/U-i / / T /\l A nv l Jl-> v " J" 

R(nv J )[l — n(nv J )\- 



A 



S(nv J — > v J ) 

"(3) 

where R(nv'J') is the excitation rate of the excited elec- 
tronic level (Equation 1), w(nv'J') is the efficiency for 
predissociation in the excited electronic state, and the 
ratio of transition probabilitie s (A nv j; Abgrall et al. 
1993a, b) is the branching ratio (|Liu fc Dalgarno| [l996). 

These models include absorption out of upper vibra- 
tional states (v > 0), and a first-order correction for self- 
absorption by H 2 at wavelengths shorter than 1120 A. 
These models also take into account electronic transi- 
tions to the B , B , D, and D states, although their 
relative contribution to the resultant spectrum is small. 
Finally, transitions to predissociating states and vibra- 
tional states that result in dissociation (v > 14, the 
vibrational continuum) are included. 

In addition to the incident radiation held, the model in- 
puts are the ground state thermal temperature (T(H 2 j), 
the column density (N(H2)), the electron energy dis- 
tribution (E e ), the Doppler b- value (in km s _1 , con- 
tributing to the line shape), and a correction for at- 
tenuation by atomic hydrogen on the incident stel- 
lar spectrum. There is evidence that the H2 tem- 
perature i n the atmosphere of HD209458b is as high 
as 2350 K ()Lecavelier P es Etan gs~et al.l 120081 ). however, 
the column density of the H 2 in the atmosphere of 
HD209458b is unknown, thus we use the Jovian system 
as a starting point for the creation of a synthetic au- 
roral/dayglow spectrum. We created a grid of models 
for several values of T{H-2), N{H2j, and E e . For the 
purposes of line identification and H2 detection, we em- 
ployed a model with T(iJ 2 ) = 1000 K (as there may 
be multiple interpretations for the data described in 
Lecavelier et al. 2008, we choose a T(H 2 ) in the range 
of the assumed effective temperature for a model exo- 
planet H 2 layer, 1000 - 1500 K; Yelle 2004; Murray-Clay 
et al . 2009), N(H 9 ) = 1 x 10 20 cm 2 (IWolven et al.l 
Il99l IWolven fc Feldmanl Il998h. and E e = 100 
eV (|Liu fc Dalgarnol 119961 ). The b- value is simply 
(2kT(H2)/m) 1 / 2 , approximately 3 km s _1 for the model 
spectrum described here, however the output spectrum is 
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Fig. 2. — Model H2 fluorescen ce spectra, bas e d on a model of the 
Jovian atmosphere presented in IWolven et al. (1997). This model 
was computed with N(H 2 ) = 10 20 cm -2 , T(H 2 ) = 1000 K, and 
E e = 100 eV. The model is described in detail in §3.1. Here we 
separate the total predicted model (blue) into components from 
electron impact-induced fluorescence (dark magenta) and fluores- 
cence excited by stellar emission lines (dark cyan), notably Lyct, 
Ly/3, C III, and O VI. The model is rebinned to 0.5 A resolution 
for display purposes. 



largely insensitive to the actual choice of b- value for val- 
ues of a few km s -1 . In the general case presented here, 
we neglect absorption by atmospheric neutral hydrogen, 
however this is most likely not proper for the HD209458b 
atmosphere, and we will return to this point in §4.2.1. 

We note that while these values most likely do not 
represent the exact conditions in the HD209458b atmo- 
sphere, the qualitative features of the output emission 
spectrum are relatively insensitive to the specifics for 
N(H 2 ) < 10 21 cm 2 and T(H 2 ) > 300 K. Very large col- 
umn densities produce enhanced emission at long wave- 
lengths (A > 1350 A) as the 912 - 1110 A flux gets self- 
absorbed and redistributed. H2 at high temperatures 
(T(i/ 2 ) ^ 2500 K) is characterized by strong band emis- 
sion in the 1200 - 1280 A range because the broadly 
distributed ground state population allows enhanced ab- 
sorption into the Werner bands which dominate the emis- 
sion spectrum following the fluorescent cascade. 

Synthetic spectra are presented in Figure 2, where we 
show a decomposition of the predicted H 2 spectrum into 
electron impact and ultraviolet line fluorescence compo- 
nents. There are several notable features of this spectrum 
in the COS (G130M + G160M) bandpass shown here: 

• Two prominent Lyman band systems roughly cen- 
tered on 1164 and 1173 A 

• Strong Werner band emission from 1200 - 1280 A 

• The "fluorescent H 2 desert" from 1300 - 1350 A 

• Several strong photon-excited emission lines in the 
1350 - 1550 region (mainly excited by stellar Lya 
and Ly/3 photons) 

• Dissociation continuum emission observed between 
1350 - 1630 A resulting from transitions to the 
unbound vibrational continuum (v > 14) of the 
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TABLE 2 

HD209458 ultraviolet emission line catalog, line 

STRENGTHS, AND LINE WIDTHS a . 





A resr 


A i 


T.irip T^lllY 


FWHM 




( A ) 




(10 — ^ ergs cm ^ s ^) 


(Km s ) 


C III b 


1176 


1174.94 


2 1Q + 31 


144.5 


C III 


1176 


1175.56 




106.7 


C III 


1176 


1176.20 


0.77 ± 0.44 


48.6 


N I 


1199.97 


1199.84 


1 00 + 30 


136.3 


Si III 


1206 50 


1206.34 


1 2 89 + 35 


73.0 


H I c ffi 


1215.67 


1215.49 


383 ± 16.71 


209.3 


O V 


1218.34 


1218.21 


3 97 + 1 5 


115.3 


N V 


1238.82 


1238.62 


1 08 + 58 


74.9 


IN V 


1Z4Z.OU 




n oft -U n O/l 
U..Z0 it U.U4 




N V d 


i O A O Oil 

1242.80 


1242.94 


1.30 ± 1.07 


134.8 


S II 


1259.52 


1259.08 


1.68 ± 0.54 


170.5 


Si II 


-1 O C A t~7 A 

1264.74 


1264.71 


2.11 ± 0.11 


117.4 


C I 


"1 O T A n O 

1274.98 


"1 n T OA 

1274.80 


0.16 ± 0.049 


28.1 


Si III 


1298.95 


1298.85 


0.63 ± 0.09 


127.8 


O I© 


1302.17 


1301.89 


7.38 ± 0.62 


70.4 


O Iffi 


1304.86 


1304.63 


11.79 ± 0.17 


133.0 


O I© 


1306.03 


1305.80 


7.42 ± 0.49 


66.8 


Si II 


1309.28 


1309.12 


o "i r" 1 n Aa 

3.15 ± 0.48 


209.6 


N I 


1 ni a nn 

1319.00 


1318.71 


0.28 ± 0.06 


21.5 


C II 


1334.53 


"1 O O /I o o 

1334.32 


o n a 1 a 

8.04 ± 0.25 


72.4 


C II 


1335.71 


1335.50 


16.06 ± 0.14 


70.2 


CI I 


1351.66 


1351.45 


"i nn 1 n on 

1.20 ± 0.30 


70.6 


O Iffi 


1355.60 


1355.39 


"i o o 1 n nn 

1.38 ± 0.09 


48.7 


Si IV 


1393.76 


1393.58 


9.47 ± 0.40 


70.1 


Si IV 


1402.77 


1402.61 


5.11 ± 0.49 


70.8 


Fe II 


1418.53 


1418.65 


0.42 ± 0.05 


41.7 


S I 


1425.03 


1424.89 


0.63 ± 0.05 


29.4 


S I 


i A o o n o 

1433.28 


1433.13 


n a ~i 1 n n o 

0.41 ± 0.08 


37.4 


C I 


1470.09 


1469.97 


0.36 ± 0.08 


36.0 


S I 


1472.97 


1472.82 


1.56 ± 0.33 


42.1 


S I 


1474.00 


1473.89 


1.33 ± 0.33 


34.6 


S I 


1481.70 


1481.54 


0.74 ± 0.22 


38.1 


S I 


1483.23 


"1 A O O AA 

1483.00 


n o 1 I f\ f\ a 

0.31 ± 0.04 


35.8 


S II e 


1519.91 


1519.78 


1.58 ± 0.87 


121.4 


Si II 


1526.71 


1526.48 


I 1 n I n n to 

1.19 ± 0.072 


23.1 


C III e 


1531.84 


1531.43 


0.75 ± 0.13 


39.1 


Si II 


1533.43 


1533.27 


2.31 ± 0.14 


63.7 


C IV 


1548.19 


1548.05 


i f\ r A 1 n o o 

10.54 ± 0.38 


72.1 


C IV 


1550.77 


1550.58 


t n*~7 1 n n/^ 

7.07 ± 0.96 


108.5 


Fe II 


1559.08 


1558.93 


1.98 ± 0.29 


61.1 


Fe II 


1563.79 


1563.61 


1.38 ± 0.74 


26.8 


Fe II 


1566.82 


1566.63 


0.73 ± 0.67 


16.0 


Fe II 


1570.24 


1570.06 


n C A 1 C\ A c\ 

0.64 ± 0.40 


26.6 


Si I+Fe II 


1573.87 


1573.71 


1.59 ± 0.10 


59.6 


Si I+Fe II 


1574.82 


1574.77 


1.67 ± 0.10 


41.8 


Fe II 


1580.63 


1580.38 


1.22 ± 0.21 


42.4 


Fe II 


1584.95 


1584.59 


1.06 ± 0.13 


43.2 


Fe II 


1588.29 


1588.09 


1.23 ± 0.13 


33.9 


Fe II 


1612.80 


1612.58 


1.61 ± 0.38 


44.9 


Fe II 


1618.47 


1618.29 


0.74 ± 0.10 


58.6 


Fe II 


1621.68 


1621.63 


0.62 ± 0.16 


38.2 


Fe II 


1623.09 


1622.88 


0.65 ± 0.06 


22.9 


Fe II 


1625.52 


1625.35 


1.47 ± 0.23 


35.8 


Fe II 


1632.67 


1632.51 


1.66 ± 0.55 


68.2 


Fel 


1637.40 


1637.26 


2.03 ± 0.17 


39.9 


He II 


1640.40 


1640.15 


2.51 ± 0.38 


63.2 


Fe II 


1643.58 


1643.38 


3.19 ± 0.26 


91.8 



a Line list created from an exposure time weighted 
coaddition of the quadrature and occultation spectra. 
b Blending of the C III multiplet led to three distinct 
fit components. 

c Lines labeled © are contaminated by geocoronal 
emission. 

d Incomplete detector pulse-height screening by the 
CALCOS pipeline leads to a spurious feature coinci- 
dent with the N V A1242 A line. 
e Tentative identification, previously undetected. 
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ground electronic state (Stecher & Williams 1967; 
Abgrall et al. 1997; see also France et al. (2005) for 
additional discussion of this feature in the context 
of photo-excited gas) 

• The strongest individual feature in this H2 spec- 
trum is the B - X (6 - 13) P(l) emission at 
1607.50 A which is excited by coincidence of the 
ground state transition B — X (6 — 0) P(l), 
1025.93 A, with stellar Ly/3 

3.2. Stellar Line Inventory 

The planetary emission signal from the HD209458 sys- 
tem is expected to be weak compared to the stronger stel- 
lar features present in the COS far-UV bandpass. As the 
far-UV spectrum of H2 can have hundreds of relatively 
strong emission lines that depend upon the specifics of 
the molecular excitation, it is critical to have a thorough 
catalog of the stellar emission features in the 1150 - 
1650 A H 2 range covered by COS. In addition, while 
high-lying ions originating in the chromosphere or tran- 
sition region of HD209458 (such as Si IV and C IV) will 
most likely not be found in an exoplanetary environment, 
it is possible that intermediate ionization species (C II, 
O I, S I, Si II; Vidal-Madjar et al. 2004) could exist in the 
hot thermosphere of a short-period exoplanet (|Yelldl2004l : 
iGarcia Munozl [20071 ) . Therefore, we compared detailed 
studies of the far-UV spectra of a Cen and Arcturus 
with our COS o bservations of HD209458 to catalog th e 
stellar features (|Pagano et al.ll2004t Iffinkle et all 12005( 1 . 
Figure 3 displays representative spectral windows from 
the G130M and G160M observations, with the H2 emis- 
sion model shown offset in blue. In Table 2, we list the 
observed emission lines, including their observed wave- 
length, integrated brightness, and line width. We note 
that stellar variability calls into question the validity of 
line fluxes measured from data acquired over the course 
of ~ 1 month. While we did find relatively good agree- 
ment between the peak line strengths at the quadrature 
and occultation pointings (± « 20%), nothing in the fol- 
lowing analysis depends on the absolute calibration of 
the stellar emission line fluxes. The line fluxes quoted in 
Table 2 would best be thought of in a relative sense. 

3.3. Quadrature Observation Cross- Correlation 

Even in the far-UV bandpass, where continuum emis- 
sion from solar type stars is weak (and consistent with 
the background at A < 1400 A), stellar contamination 
must be addressed. Steps must be taken in data analysis 
to maximize the potential emission from the planet while 
minimizing the contribution from stellar line and contin- 
uum emission. We take advantage of the orbital velocity 
of HD209458b, v or b = 144 km s assuming a circular or- 
bit, and coadd the quadrature observations as a function 
of velocity. Additionally, we subtract the stellar emission 
observed at the occultation (Phase 0.5) position in order 
to remove additional confusion from weak stellar lines. 
Stellar variability is accounted for with a scaling based 
on the flux in high- ionization chromospheric/transition 
region lines (e.g. Si IV A1400 and C IV A1550 A) not 
expected to be present in the planetary atmosphere. We 
found that HD209458 was relatively stable across all of 




1220 



1570 1580 1590 1600 1610 1620 
Wavelength (A) 

Fig. 3. — Close-up views of individual spectral windows observed 
with COS G130M (top) and G160M (bottom). Prominent stellar 
and geocoronal features are labeled. The model H2 spectrum de- 
scribed in 53.1 is shown offset in blue. 



the observations used here (Section 3.1), with peak vari- 
ability of AF(C IV) < 5 x 10~ 16 ergs cm" 2 s" 1 A" 1 , 
which is « 19% of the stellar flux at this wavelength. 

This process not only increases the signal from the 
planet, but serves as a cross-correlation function that al- 
lows one to assess the validity of an emission signal. We 
would expect any signal that originates on HD209458b 
to be at roughly the orbital velocity of the planet. Math- 
ematically, we define a star-subtracted quadrature spec- 
trum (D(v or b,X)) which is a function of both wavelength 
and the orbital velocity of the planet, as 



D(vorb,\) = (F Q1 {X)+F Q2 (X+dX))-{F ccW- 
with 



dX 



2v n 



-A„ 



-F 0c c(X+dX)) 
(4) 

(5) 



and A c = 1400 A, where F Q1 , F Q2 , 
fluxed spectra (in units of ergs cm~ 2 



and Focc are the 
s _1 A -1 ) observed 

at the first quadrature (Phase 0.25), second quadrature 
(Phase 0.75), and occultation (Phase 0.5), respectively. 
We create this planetary spectrum as a function of or- 
bital velocity for the range 50 < v or b < 250 km s _1 , in 
1 km s _1 bins. In order to suppress noise spikes intro- 
duced during the subtraction, we rebin D(v or b,X) to 0.1 A 
pixel -1 (« 20 km s -1 ), though numerous binnings were 
tested in order to rule out spurious features introduced 
by a particular choice of pixel scale. The correlation sig- 
nal in this treatment is the peak flux of D(v or b 7 X) in a 
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given velocity bin and spectral range, 

C(v orb ) = MAX(D(v or b, AA)) (6) 

where AA is the wavelength interval around the emission 
line of interest, typically ~ 10 A. 

Several spectral regions were identified to isolate spe- 
cific emission lines of interest that fell into one of the fol- 
lowing categories: 1) emission line-free, a region contain- 
ing stellar continuum or detector background, 2) emis- 
sion line regions containing only stellar features, 3) emis- 
sion from low-ionization ions and neutrals which could 
have components from both the star and the planet, and 
4) regions that should contain a wealth of planetary H2 
features. A full list of the spectral windows investigated 
is provided in Table 3. We then look for peaks in the 
correlation signal, C{v or b), as a function of orbital ve- 
locity. Figure 4 displays the results for six of these re- 
gions. We also measure the wavelength of peak emission 
in order to verify that the signal is associated with a 
specific line and not a chance noise peak. As an illustra- 
tive comparison, we describe plots (a) and (b) in Figure 
4. These are regions dominated by stellar continuum 
(1500 - 1510 A) and stellar emission lines (C IV doublet, 
1545 - 1555 A) respectively. The left panels show the 
correlation peaks (C(v or b), in units of 10 -15 ergs cm~ 2 
s _1 A -1 ) as a function of orbital velocity. Mean fluxes 
and standard deviations are computed from an average 
of two velocity regions distinct from the circular velocity 
(60 < v orb < 100 and 200 < v orb < 240 km s" 1 ) and 
displayed in Figure 4. While the C(v or b) plots show a 
mostly random distribution, the right panels show that 
while the peak wavelength of the continuum region is 
approximately randomly distributed (a) , residuals in the 
fluxes in the C IV 1550.77 A are responsible for most of 
the peaks in the 1545 - 1555 A region (b). 

Figure 2 (c - f) shows four additional spectral regions 
that highl ight possible emissi o ns fro m the HD209458b 
system. IVida l-Madi ar et all (|2004l ) observe possible 
absorption in the C II 1335 A multiplet, arising in 
the extended envelope of the planet, therefore this 
could b e a viable emitting species in a hot t hermo- 
sphere (jGarci'a Muiiozl 120071: iLinskv et al.l 12010ft . O I 
1356 A is considerably less contaminated by geocoro- 
nal emission than the O I 1304 A multiplet, and could 
be present in the exoplanet. In fact, the O I 1356 A 
emission is observed to be brighter than O I 1304 A 
in the Jovian system, where O I 1356 A is produced 
throu gh the dissociation of O2 and H2O by electron im- 
pact |Hall et all 1 19981: iNoren et al.l 120011) . In panels (c) 
and (d), the peak wavelength, Ac, displays a linear in- 
crease with orbital velocity, that we attribute to a broad 
geocoronal component that is being shifted to longer 
wavelengths in this treatment. In the bottom panels 
(e and f), we display the correlation peaks for two re- 
gions with an abundance of strong H2 featu res, centered 
around 1575 and 1608 A (|Gustin et al.ll2009ft . Nothing is 
seen in the 1608 A band, but there are several correlation 
peaks present in the 1575 A region within ± 20 km s _1 
of the planet's nominal circular velocity. The strongest 
two include a broad feature, approximately centered on 
the circular velocity of the planet, and the second is lo- 
cated at 158 - 160 km s" 1 (w +15 km s _1 in the rest 



frame of the planet). We note that HD209458b is nearly 
circularized (e = 0.014 ± 0.009, consistent with a circu- 
lar orbit; Laughlin et al. 2005), and the amplitude of 
variation in the orbital velocity is of order a few km s . 
The significance of these results and the possibility of an 
H2 detection are discussed in §4. 

4. RESULTS AND DISCUSSION 
4.1. Emission Feature at 1582 A 

An analysis of C(v or b) for each of the spectral sub- 
regions identified in Table 3 indicates two statistically 
significant features that are consistent with emission 
from HD209458b. There is a strong peak in the cross- 
correlation spectrum at v orb = 158 - 160 km s _1 in the 
1575 A H2 band. This 1582 A feature is detected at the 
3.8 a level (> 99.9% confidence) at 1581.95 ± 0.10 A 
(Figure 5). The bottom panel of Figure 5 shows 
the quadrature summed, stellar subtracted spectrum 
(D(v or b,X); Equation 4) for v or b = 159 km s _1 . This fea- 
ture is not obviously detected in either quadrature spec- 
trum individually, and does not correspond to a known 
stellar feature. In Figure 3, we showed that there are 
several Fe II lin es in this region. T he -ffST-STIS spec- 
trum of a Cen ([Hinkle et al.l I2005T ) contains weak emis- 
sion from Fe II A 1581.27 A. While we cannot explicitly 
rule out this line as the explanation for the 1582 A fea- 
ture seen in HD209458b, we consider it unlikely that only 
this Fe II line is present when dozens of stronger Fe II 
lines are not observed. 

We calculated the probability of randomly finding a 
3.8 a event in the 166 A sampled in the 14 intervals 
(1660 total samples at 0.1 A binning) listed in Table 3. 
Assuming Gaussian statistics, we find that one should 
expect 0.4 such events in the entire range, for each ve- 
locity interval. This analysis indicates two (at v or b = 158 
and 159 km s^ 1 ) 3.8 a events, from which we estimate 
the possibility of a false positive at ~ 16%. It should be 
emphasized that this is the probability of finding a spuri- 
ous signal anywhere in the 166 A sampled. In reality, the 
false positive probability is much lower because there are 
only a small range of velocities that could be plausibly 
consistent with an ex oplanetary origin. Given the ~ 42 
km s" 1 escape speed (jLinskv et al.ll2010l ). a conservative 
range of possible atmpsheric emission velocities is ± 60 
km s _1 , or ± ~ 0.6 A, using the binning described above. 
This reduces the possibility of a random detection in each 
velocity interval by a factor of 20, leading to a realistic 
false positive probability of ~ 0.04%. Thus, we conclude 
that this feature, at +15 km s _1 in the rest frame of the 
planet, is most likely a real signal in the data. 

There are several strong H2 lines in this region, the 
brightest being the Lyman band transitions (B - X) 
of (2 - 8) P(12), (7 - 13) P(3), (6 - 12) P(l), (2 - 
9) P(4), and (7 - 13) P(4), X rest = [1580.70, 1580.74, 
1581.11, 1581.11, 1581.39 A], respectively. The oscillator 
strengths for each of the transitions are larger than 0.063, 
except for (6 - 12) P(l) which has / = 0.027. However, 
(6 - 12) P(l) has a pumping transition coincident with 
Ly/3, as discussed in the following section. The (7 - 13) 
lines are excited by electron impact, and are not seen else- 
where in the spectrum. In the bottom panel of Figure 5, 
we overplot the UV photon fluorescence-only H2 model 
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Fig. 4. — Plots of the correlation peak, C{v or b) (defined in §3.3) for orbital velocities from 50 - 250 km s _1 . Several representative 
spectral regions are displayed: (a) stellar continuum only, (b) stellar emission lines, (c) and (d) stellar atomic features with potential 
contributions from HD209458b, and (e) and (f) H2 band regions. For each v or h, the plots on the right show the wavelength corresponding 
to the maximum of D(?) or j,,A). The dashed gray vertical line represents the rest frame of the exoplanet for a circular velocity of 144 km 
=-1 



Auroral/Dayglow Observations of HD209458b 



TABLE 3 

QlJADTRATURE OBSERVATION EMISSION CORRELATIONS. 



Orbital Velocity (km s ) Orbital Velocity (km s ) 

120 140 160 180 120 140 160 180 



Species 


Origin a 


AA 


Peak Flux Limit b 






(A) 




Background 




1145 - 1155 


1.34 


C III 


ST 


1172 - 1180 


1.30 


H I 


ST, PL, ® 


1212 - 1222 


1.14 


Background 




1280 - 1290 


0.48 


Ol 


ST, PL, ® 


1300 - 1310 


2.56 


C II 


ST, PL 


1330 - 1340 


2.92 


O I 


ST, PL 


1350 - 1360 


1.21 


Si IV 


ST 


1390 - 1406 


2.32 


S I 


ST, PL 


1470 - 1478 


0.99 


Continuum 


ST 


1500 - 1510 


1.11 


Si II 


ST, PL 


1524 - 1538 


1.31 


C IV 


ST 


1545 - 1555 


1.79 


H 2 +Si I+Fe II 


PL, ST 


1565 - 1585 


2.05 c 


H 2 +Fe II 


PL, ST 


1600 - 1620 


3.26 




a Proposed physical origin in the HD209458 observations: 

ST = Stellar, PL = Planet, ffi = Geocoronal airglow. 

b in units of (10 — 16 ergs cm -2 s^ 1 A -1 ). 1 <r limit on the orbital 

velocity cross-correlation, as described in §3.3, computed over 60 

< v orb < 100 and 200 < v orb < 240 kms' 1 . 

c Tentative H2 detection at 3.8 a above this limit, discussed in §4. 

in blue, shifted to v orb = 159 km s" 1 . The (2 - 8) P(12) 
line, which is pumped by coincidence with C III 1175, is 
not seen. Both of the 1581.11 A transitions are shifted 
to 1581.95 A at the velocity of the strong correlation 
peak. The other features are displaced by |AA| > 0.25 A, 
larger than the COS resolution element, and the 0.1 A 
binning applied to the D(v or b,X) spectrum. As noted in 
§3.3, D(v or b,X) was computed for several binnings to test 
for spurious features. As would be expected for an un- 
resolved molecular line, this feature is not detected for 
rebinning > 1 A. We conclude that B - X (6 - 12) P(l) 
A1581.ll A and B - X (2 - 9) P(4) A1581.ll A are the 
most likely candidates for the possibly detected feature. 
In Section 4.2 and 4.2.1, we will argue that this feature, 
if real, cannot be Ly/3-pumped (6 - 12) P(l), and we 
describe a possible excitation mechanism for the (2 - 9) 
P(4) line. 

4.2. H 2 Emission from HD209458b 

Strong auroral and/or dayglow H2 emission from 
HD209458b is not observed. Figure 3 is representative 
of the results from comparisons between our COS obser- 
vations and the models described in 3.1. Upper limits 
on the line strengths of the strongest emission lines pre- 
dicted by our model are presented in Table 4. As dis- 
cussed above, the two candidate emission lines for the 
observed feature in the D(v or b,X) spectrum are (6 - 12) 
P(l) and (2 - 9) P(4). This feature peaks at +15 km s" 1 
in the rest frame of the planet, with a second broader 
maximum approximately at rest in the planet frame (Fig- 
ure 5). The (2 - 9) P(4) line is intrinsically stronger 
than the (6 - 12) P(l), however a lack of clear excitation 
mechanism complicates this interpretation of the result. 
There is no obvious pumping transition ((2 - v) R(2) or 
P(4)) exactly coincident with an observed stellar emis- 
sion line. Possibilities exist if there is a velocity offset 
between the "H2 layer" of the planet (1.0 < r < 1.1 Rp, 
where Rp is the radius of HD209458b; Murray-Clay et 
al. 2009). For example, the (2 - 0) band is near strong 
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Fig. 5. — Same as Figure 4, but showing the cross-correlation 
spectra in greater detail. One observes three peaks found at greater 
than 2 a significance within ± 20 km s — 1 of the planet rest frame 
( v orb = 144 km s _1 ), and a peak at 158 - 160 km s~ 1 detected at 
3.8 a confidence. The difference spectrum, D{v orb ,X), at v orb = 159 
km s — 1 is shown in the lower panel, with the UV photon fluores- 
cence model shown overplotted. Both B — X (2 — 9) P(4) and 
B - X (6 - 12) P(l) emit at A rest = 1581.11 A, however, the 
absence of additional lines from the Ly/3 progression (described in 
§4.2.1) indicate that the observed line is B - X (2-9) P(4), which 
may be excited by stellar O I emission (§4.2). 



N II and He II tra nsitions (A 1083 - 10 86 A) observed 
in solar-type stars (jRedfield et al.l 1200 2). We analyzed 
FUSE observations of HD209458 in order to constrain 
the stellar spectrum in this wavelength region, but the 
data were not of sufficient quality to derive meaning- 
ful constraints. Another possibility is coincidence be- 
tween (2 - 4) P(4) 1301.70 A and O I 1302. For 6- 
function emission and absorption profiles, these lines 
would require a velocity offset of 108 km s _1 to over- 
lap, however in practice these lines will have an appre- 
ciab le velocity width that may prov ide sufficient over- 
lap (Ben-Jaffel & Son a Hosseinil [20091 ) . 

Exploring the possibility of O I excitation further, 
we examined the relative line fluxes of stellar Lya and 
O I. The data set for HD209458 is complicated by ob- 
servations with different instruments, spectral resolu- 
tions, and the attenuation of both lines by 47 pc of 
ISM. A somewhat more straightforward determination 
can be ma de using the high-re ssolution STIS spectrum 
of a Cen (Pa gano et al.l [20041 ) . Interstellar absorption 
still complicates the comparison, but we estimate an 
F(Lya)/F(0 I) ratio of - 20 - 40 for HD209458. When 
considering the possibility of H2 excited by stellar O I, 
it is important to note that appreciable occupation of 
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TABLE 4 

Limits on H2 emission from the atmosphere of HD209458b. 





Line ID b 


(A) 


Excitation 


Line Flux 
(10 ergs cm s ) 


Notes c 


c 


- X (1 - 4) P(5) 


1171.96 


e — Impact 


< 0.76 




c 


-X (2-5) Q(l) 


1175.83 


e~ Impact 


< 6.22 


Blend with C III 




-A (4-8) Q(l) 


1 nan r o 
1239.0O 


e — Impact 


< 0.44 




B 


- X (1 - 3) R(5) 


1265.67 


UV photon, Lya 


< 0.71 




C 


- X (7 - 12) Q(3) 


1279.10 


e — Impact 


< 0.80 




B 


- X (6 - 7) P(l) 


1365.66 


UV photon, Ly/3 


< 2.01 




B 


- X (6 - 9) P(l) 


1461.97 


UV photon, Ly/3 


< 3.36 


Stellar continuum at A > 1420 A 


B 


- X (1 - 8) R(3) 


1547.34 


UV photon, Lya 


< 5.83 




B 


— X (7 — 13) P(3) 


1580.74 


c~ Impact 


< 10.14 


Blend with Fe II 


B 


- X (2-9) P(4) 


1581.11 


UV photon 


11.72±3.08 


Blend with Fe II 


B 


- X (6 - 12) P(l) 


1581.11 


UV photon, Ly,3 




Blend with Fe II 


B 


- X (6 - 13) P(l) 


1607.50 


UV photon, Ly,3 


< 10.54 




B 


- X (5 - 12) P(3) 


1613.18 


e — Impact + 

TTV Plin+nr, n vt 


< 12.80 





UV Photon, O VI 



the v = 4 level of the ground electronic state requires an 
energetic population of H2, but this may be reasonable 
to expect in the atmosphere of a hot Jupiter. The at- 
mospheric temperature in the HD209 458b atmosphere i s 
expected to rise sharply at r > Rp (Sin g et al.ll2008bh . 
and this could increase the rovibrational temperature of 
the molecules in the H2 layer with height, up to the col- 
lisional dissociation threshold for H2 (T « 4000 K; Shull 
& Beckwith, 1982). Additionally, there is evidence that 
the vibrational temperature of H2 in the Jovian atmo- 
sphe re can exceed the kinetic temperature by as much as 
50% (fBarthclc mv et aL| [2005) and electronic excitation of 
H2 from highly excited levels of the ground state is ob- 
served in the warm H2 environment around the classical 
T Tauri star TW Hya (jHerczeg et al.ll2002D . 

4.2.1. Constraints on Surface HI Lyman Series Strength 
and Magnetic Field 

Given the uncertainties in the excitation route for the 
(2 - 9) P(4) emission line described above, it would 
be tempting to assign the 1582 A emission to (6 - 12) 
P(l) pumped by stellar Ly/3 through the (6 - 0) P(l) 
1025.93 A transition. This scenario does not hold as the 
(6 -12) band of the Ly/3 pumped progression is among 

the weaker branches of (6 - v ). There are several other 
lines excited by stellar Ly/3 that should be > 2 times 
brighter (and thus detectable), including the (6 - 13) 
P(l) A1607.50 A line that should be ~ 4 times brighter 
than (6 - 12) P(l). As none of these lines are observed, 
we are forced to rule out Ly/3 pumping for the observed 
feature. 

The H I Lyman series lines are the strongest far-UV 
spectral features in a solar-type star, and assuming that 
the 1582 A feature is produced by (2 - 9) P(4), this 
raises the question of why H2 lines excited by L yman 
photo ns are not observed. iVidal-Madiar et all (2003, 
120041 ) describe the detection of an extended H I atmo- 
sphere surrounding the planet, and we hypothesize that 
this warm cloud of neutral hydrogen is primarily respon- 
sible for the suppression of Lyman-pumped H2 emission 
from HD209458b. This H I layer will scatter the in- 
cident Lyman series lines, effectively diffusing the pho- 
tons away from the line core where the H2 coincidence 
exists. By comparing model line-strength ratios with 



the observed F(A rest 1581.11 A)/F(X rest 1547.34 A) and 

rest 

1581.11 k)/F(X 

rest 

1607.50 A) ratios, we can con- 
strain the relative optical depths seen by stellar Lya and 
Ly/3 photons before they reach the H2 layer. We take 
the 1 a limits for the "measured" fluxes of the 1547.34 
and 1607.50 A lines, and the model line ratios are « 0.47 
and 0.25, respectively. Assuming that attenuation of the 
pumping transition for (2 - 9) P(4) is negligible, we es- 
timate that optical depths at the (1 - 2) R(6) 1215.73 A 
and (6 - 0) P(l) 1025.93 A transitions are 8.1 and 15.2, 
respectively. We created a simulated H I profile to deter- 
mine the column density required to attenuate the stellar 
Lyman lines, finding N(HI) > 3 x 10 20 cm~ 2 on the line 
of sight from HD209458 to the H 2 layer on HD209458b. 

The far-UV spectra of Jupiter and S aturn are dom- 
inate d by electron impact excitation dFeldman et al.l 
119931: iLiu fe Dalgarnol Il996t iGustin et all I2009D . and 
in addition to understanding the lack of H2 emission 
pumped by Lyman series photons, the apparent lack 
of electron impact H2 must be addressed. The ener- 
getic electron spectrum is driven by the interaction of 
the planetary magnetic field with the ambient atmo- 
sphere (Broadfoot et al. 1979, and see Bhardwaj & 
Gladstone 2000 for a comprehensive review of the pro- 
duction of H2 emissions in the Jovian aurorae). Thus, 
assuming that the atmosphere of HD209458b has an H2- 
dominated composition similar to that of Jupiter, the 
paucity of electron impact emission can be interpreted 
as a reflection of a relatively weak surface magnetic field 
on HD209458b. This interpretation is supported by 
theoretical studies of the magnetic environment of hot 
Jupiters. The magnetic moment (A4) is known to scale 
with the core d ensity of the planet, t he rotation rate, 
and core radius (|GrieBmeier et aLll2004D . A short-period 
exoplanet such as HD209458b is expected to be tidally 
locked due to gravitational dissipation, which results in a 
dramatically lower rotation rate compared to a gas giant 
such as Jupiter. Combining the rotation rate, the den- 
sity of the planet inferred from transit observations, and 
models of the interior structure, HD209458b is expected 
to have an intrinsic magnetic moment (in the tidall y 
locked limit) of M < 0.01 M j (|Griefimeier et al.H2004D . 
When factoring in the interaction with the stellar mag- 
netic field, this value is of order 0.1 A4j, consistent with 
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the su rface magnetic field calculations of San chez-Lavegal 
(|2004D . B < 0.1 Bj (where Bj is the Jovian polar 
magnetic field). Interestingly, recent results from radio- 
frequency o bservations of the well-studied ho t Jupiter 
HD189733b (jLecavelier Pes Etangs et all 120091) support 
the hypothesis of a weak exoplanet magnetic field. An- 
other possible source of energy for atmospheric elec- 
trons is the interplanetary magnetic field (-Bjp), but 
this is estimated to b e small for the HD 209458 system 
B IP < 5 x 10~ 3 Bj (jErkaev et all [2001 . We conclude 
that the feature observed at 1582 A in HD209458b is 
consistent with no detectable lines from electron impact 
H-2. The slow rotation and low density of the planet cre- 
ates a magnetic field that is incapable of producing H2 
emission with intensities observed in the Jovian system. 

4.2.2. Alternative Excitation Routes and H2 Radiative 
Transfer Effects 

We have shown that the primary excitation routes for 
populating the Lyman and Werner bands in gas giants in 
our solar system may be suppressed in the atmosphere of 
HD209458b. It is worth noting that the higher lying elec- 
tronic states of H2 can be populated via electron impact 
and absorption of extreme UV (EUV; A < 911 A) pho- 
tons. Electron impact excitation of the E and F states 
of H2 can populate the Lyman bands via a radiative cas- 
cade process (E, F 1 T,+ — > B 1 S+ + hv), contributing as 
much as ~ 25% of the tota l flux from the B - X band 
emission ( Aie llo et al.l lT982'). We note that this process 
concentrates the emitted photons in the 1350 - 1500 A 
region, and does not favor the putative B - X (2 - 9) 
P(4) emission feature. Coupling this with the arguments 
against electron impact excitation noted above, we dis- 
miss this as a possible energy source. What we cannot 
rule out is the possibility that one of the higher-lying 
electronic states is photo-excited by EUV stellar line pho- 
tons. This cascade could preferentially populate a spe- 
cific level of the B state, causing the resultant emission 
spectrum to differ from that expected from direct exci- 
tation of the Lyman and Werner bands. However, the 
upper H I atmosphere of HD209458b should effectively 
attenuate any stellar EUV emission capable of exciting 
H2, making this route unlikely. 

One possibility that could explain the lack of H2 fea- 
tures seen in the spectrum of HD209458b is the coherent 
scattering of pumping photons absorbed in the wings of 
the H2 line p rofiles. This par t ial frequency redistribu- 
tion process (jLupu et al.l l2009t iBarthelemv et al.l 120041 ) 
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requires a more sophisticated radiative transfer treat- 
ment than the synthetic H2 models presented in this 
work, where all photons absorbed in a given transition 
are redistributed to the line core prior to emission. This 
partial redistribution process, which has been observed 
in the spectra of Jupiter and Saturn, acts to suppress 
the line-to-continuum ratio and decrease the contribu- 
tion from subordinate transitions in a given progression. 
Additional radiative transfer investigations into EUV ex- 
citation of H2 and frequency redistribution would be very 
useful for the interpretation of future experiments char- 
acterizing the far-UV spectra of extrasolar giant planets. 

5. SUMMARY 

We have presented new HST-COS observations aimed 
at characterizing aurorae on the extrasolar planet 
HD209458b. We use a model for the auroral and day- 
glow spectra of Jupiter to predict the expected features 
and use this spectrum to guide the data analysis. These 
data reveal no unambiguously identifiable signature of 
atomic or H2 emission from the planet. We do detect a 
significant feature at 1582 A, observed at +15 (± 20) km 
s _1 in the rest frame of the planet. This feature may be 
the H 2 B - X (2 - 9) P(4) line at A rest = 1581.11 A, 
however the excitation mechanism for this line is un- 
clear. Based on this feature, we placed limits on the level 
of emission from electron impact and stellar fluorescent 
excitation. We put these limits in the context of the 
HD209458 system, deriving a limit on the atmospheric 
column of H I. The non-detection of electron impact H2 is 
consistent with predictions of a weak magnetic field due 
to tidal locking and the low density of the planet. We 
suggest that future ultraviolet studies of auroral/dayglow 
emission from extrasolar giant planets would benefit by 
targeting the nearest transiting planet systems that do 
not show evidence for extended atomic and ionic atmo- 
spheres. 
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